Non-Markovian dynamics in pulsed and continuous wave atom lasers 
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The dynamics of atom lasers with a continuous output coupler based on two-photon Raman 
transitions is investigated. With the help of the time-convolutionless projection operator technique 
the quantum master equations for pulsed and continuous wave (cw) atom lasers are derived. In the 
case of the pulsed atom laser the power of the time-convolutionless projection operator technique 
is demonstrated through comparison with the exact solution. It is shown that in an intermediate 
coupling regime where the Born-Markov approximation fails the results of this algorithm agree with 
the exact solution. To study the dynamics of a continuous wave atom laser a pump mechanism is 
included in the model. Whereas the pump mechanism is treated within the Born-Markov approx- 
imation, the output coupling leads to non-Markovian effects. The solution of the master equation 
resulting from the time-convolutionless projection operator technique exhibits strong oscillations in 
the occupation number of the Bose-Einstein condensate. These oscillations are traced back to a 
quantum interference which is due to the non-Markovian dynamics and which decays slowly in time 
as a result of the dispersion relation for massive particles. 



I. INTRODUCTION 



Nowadays it is a standard technique toproduce a Bose- 
Einstein condensate in the laboratory [|l]||]. In order to 
build a coherent source of atoms, an atom laser, a major 
achievement was the coherent extraction of atoms from 
an atomic trap. At first a pulsed atom laser was built 
HH , recently also continuous wave atom lasers have been 
realized [plpj . A short survey over the experimental sit- 
uation is given in |j,g] . 

The theoretical treatment of an atom laser is usu- 
ally based on the Born-Markov approximation [p|dl3[, 
which has been used in quantum optics with great suc- 
cess. However, as it has been shown recently by Moy 
and coworkers |14[] , this approximation fails in a realistic 
parameter regime. In this article we will outline a differ- 
ent approach, which is based on the time-convolutionless 
(TCL) projection operator technique |l5|-|l7|] to study 
non-Markovian effects resulting from the output cou- 
pling. This technique is based on a perturbative expan- 
sion in powers of the output coupling strength. 

The master equation resulting from the perturbative 
expansion has a form similar to the Born-Markov master 
equation and is also local in time. This makes the equa- 
tion of motion easy to solve. As we will see the second or- 
der perturbation theory corresponds to the Born-Markov 
approximation. By taking higher orders of the expansion 
non-Markovian effects can be studied in a systematic way. 
It is demonstrated that the perturbative expansion holds 
in an intermediate coupling regime, which corresponds to 
realistic parameters, while the Born-Markov approxima- 
tion fails for these parameters. 

This paper consists of two parts and is structured 
as follows. In section O we investigate the validity of 
the time-convolutionless projection operator technique 
(TCL) for a pulsed atom laser. The model investigated 
was discussed by Moy et al. H] . They showed that the 



Born-Markov approximation fails for realistic parame- 
ters. We demonstrate that the time-convolutionless pro- 
jection operator technique agrees with the exact solution 
for these parameters. Afterwards we apply the TCL al- 
gorithm to a simple model of a continuous wave atom 
laser in section [II. The numerical results obtained from 



a simulation using a perturbation expansion including 
4th order clearly reveal strong oscillations in the occupa- 
tion number of the Bose-Einstein condensate (BEC). As 
will be shown these oscillations can be interpreted as a 
quantum interference effect which clearly reveals depar- 
tures from the golden rule and demonstrating the non- 
Markovian dynamics of the atom laser. A short summary 
concludes the article. 



II. PULSED ATOM LASERS 



In this part we will investigate the dynamics of a pulsed 
atom laser. This means that there is an output coupling 
mechanism which transfers the atoms out of the trap but 
there is no pump mechanism which supplies the BEC 
with new atoms. When the trap is empty it is replenished 
and a new cycle starts. Therefore our starting point is 
an existing BEC inside an atomic trap and we study only 
the output coupling. 

This part is organized as follows. We briefly intro- 
duce our atom laser model and derive the exact equation 
of motion for the expectation value of t he atomic num- 
ber operator of the BEC in section II A . In section [[B 
we discuss the Born-Markov approximation. The TCL- 



algorithm is then applied to the model in section [[ C . 
Afterwards the numerical results are discussed in sec- 
A model of a continuous wave atom laser which 
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includes a pumping mechanism to compensate losses by 
the output coupling will be introduced in part III. 



A. Exact solution 

Usually, the output coupling from a BEC to a large 
reservoir is described within the Born-Markov approxi- 
mation. This description is based upon a master equa- 
tion of the form 



dt 



Psys(i) = 1m V[a] p S y S (i) 



(1) 



where the single trap mode occupied by the BEC is de- 
scribed by the creation and annihilation operators a) , a 
and 7m is the Markovian decay rate. The superoperator 
T> is of Lindblad form Eq] and is defined by 



1 



V[a] P = aPa) - -(a 1 aP + Pa) a 



(2) 



In the Markovian approach one obtains a master equation 
containing only system variables. Two basic assumptions 
underlying the Born-Markov approximation are that (i) 
if the system and the reservoir are uncorrelated at the 
beginning they remain uncorrelated for later times and 
(ii) that the evolution of the system is Markovian. The 
Markovian property means that the future evolution of 
the system depends only on the present state and is in- 
dependent of the previous history of the system. As it 
is shown in Ref. [Of the Born-Markov approximation is 
not valid in a realistic parameter regime of atom lasers. 
We want to investigate a Raman output coupler 
through state change, which was suggested by Moy et 
al. 



[14 



|21| . Two lasers are tuned to a two-photon res- 
onance to couple an initial atomic state inside the trap 
to a final atomic state outside the trap. Due to the con- 
servation of linear momentum the atoms receive momen- 
tum kicks during the Raman transitions which pushes 
the atoms out of the trap. If we assume that the lasers 
are detuned far from single photon resonances, then all 
initially empty modes of the atomic trap remain empty 
and we can neglect all modes besides the ground state 
mode occupied from the BEC. In addition, we also ig- 
nore the effects of atom-atom interactions. The resulting 
Hamiltonian H is of the form 



H — H sys + H rcs + H mt , 



H sys = Tiloq a) 



a. 
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H int = -ih{Bo) - B^a) , 



where we have defined 

f + OO 
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2M 



(4) 



The bosonic creation and annihilation operators for the 
free atomic state are denoted by b\, &&. The ground state 
trap energy is denoted by Tiujq and TiuJk is the energy of 



the free output atomic state, M is the atomic mass. The 
function n{k) describes the strength and spectral form 
of the coupling. In the specific case of a Raman output 
coupler and a harmonic trap with a Gaussian ground 
state pCipiJ of width a k in fc-space given by 



(5) 



*(*) = (2na k )- 1 / 4 exp(-fc 2 /(4^ 2 )) 



one obtains together with Eq. (3d) for the interaction 
Hamiltonian H m t 



K(k) = ^H= exp(-(fc - k f/{Aal)) . 



V^l 



(6) 



Here T 1 / 2 denotes the coupling strength of the output 
coupling and hko is the momentum transferred through 
the Raman transition. For the sake of simplicity we as- 
sume that we have two counter propagating laser beams 
with the wave vectors fe^i ~ — &L2- Therefore we have 
fco ~ 0. With the atomic dispersion relation from Eq. 
(f|) one derives the spectral density J(w) of the output 
coupling strength, 



J(w) = 



r 



,—u/a 



y/TTOUJ 

which diverges at oj = 0. The quantity a is given by 

a= 2i- 



(7) 



(8) 



With the help of the function J(us) we can deduce the 
Heisenberg equations of motion for the system operators 
a^ and a^a. Assuming that the reservoir is empty at the 
beginning, (b' k bk}(0) — 0, we obtain 



dt 



tf(t))=iLJ (a!(t)) 



dr/*(T)(ot(t-r))( 



(9a) 







Tr{p rcs B(r)i?t(0)} e ^ 
dcj J(uj)e l ^-^ T 



j t {aHt)a(t)) =- J dr f*(r){a^t) a(t - r))e~^ 

+ h.c. (9b) 

(10) 
(11) 

-. (12) 
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is the reservoir correlation function apart from a factor 

The formal solution of the integro-differential equa- 
tions (||) can be expressed in terms of inverse Laplace 
transforms. However, the numerical evaluation of the in- 
verse Laplace transforms is very difficult p3] . As is easily 
shown, the solution of Eq. (9b) can be written as 



(3a) 


The function 


(3b) 


f(r) 


(3c) 




(3d) 









(a t a)(t)=c(t)c*(t) 



(13) 



where c* (t) is a solution of Eq 
c(0) 



Pq ) with the initial value 
Hence one solves Eq. (Bah through 



vW(o) 

direct numerical integration and obtains the expectation 
value of the atom number by taking the squared absolute 
value. 



B. Born-Markov approximation 

Now we apply the Born-Markov approximation to the 
above model of an atom laser. The Born-Markov approx- 
imation consists in making the Born approximation that 
is, one assumes that system and reservoir are uncorre- 
cted initially and that they remain so for later times. 
Hence the total density operator ptot(i) can be written 



Ptot(i) = Psys(t) <8> Pres(O) , 



(14) 



where p sys (t) denotes the system density operator. In 
addition, it is assumed here that the reservoir density 
operator does not change in time p rcs {t) = /Orcs(0). The 
Markov approximation is based on the assumption that 
the time scale t sys , which describes the relaxation of the 
reduced system, is much greater than the time scale t res 
which represents a measure for the width of the reservoir 
correlation function, i. c, 



^sys ^ > £r 



(15) 



The time scale i ros is obtained from the reservoir corre- 
lation function /(r) . Because of the slow decaying shape 
of /(r) in Eq. (n2f) the atom laser shows a strong non- 
Markovian behavior. A more detailed discussion of the 
application of the Born-Markov approximation to atom 
lasers can be found in Ref. |IJ,^|. Within the Born- 
Markov approximation we get from Eqs. (pft the follow- 
ing master equation: 







-^rPsysW = -7; s MWa, p sys {t)) + j M < —xQ yap sys (t) 



dt 



1 



--p sy a(t)a^a + ap ays (t)a^ 



(16) 



Here, the Markovian Lamb shift Sm and the Markovian 
decay rate jm take the form 

/■+O0 /-+oo 

1m = / dt <j>{t) , S M = / dt V-(i) . (17) 
Jo Jo 

The functions <p(t) and tp(t) are proportional to the real 
and imaginary part of f(t), 



/(<)=oW*)+^(<)) 



(18) 



From the master equation (Hq) we easily get the expec- 
tation value of the atom number operator a^a 



(a t a)(t) = (a t a)(0)e- 7M * 



(19) 



The above results enable one to derive from inequality 
( |l5| ) an explicit condition for the Born-Markov approxi- 
mation. The system time scale t sys is given by 



1 

1m 
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(20) 



This result is easily obtained by performing the first in- 
tegral in ([l7|). As discussed in Q we define t res as the 
half- width of the integral of the real part of /(r). The 
resulting equation is solved numerically. In the consid- 
ered parameter regime one obtains i res sw OA/u>o. Hence 
the time scale condition 



t 



«y« 



3/2 
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(21) 



must be satisfied for the Born-Markov approximation to 
be valid. 

In our simulations we take similar parameters as dis- 
cussed in @. The atom mass is M « 2 X 1CT 26 kg. 
Realistic parameters for the system frequency lvq and the 
coupling strength T are w /(2tt) rs 123 s^ 1 , T w 10 5 s" 2 
P-p|, p0p^ ] . The standard deviation a^ in fc-space of the 
coupling function re(fc) is assumed to be <7k ~ 10 6 m _1 
corresponding to a wave length A ss 2/jm. 

The only parameter which is varied is the coupling 
strength T, thereby we can systematically change the ra- 
tio of the time scale condition (21). In Fig. [I] we depict 
the normalized expectation value of the atomic number 
operator obtained from the exact solution and from the 
Born-Markov approximation for V = 5 x 10 4 s~ 2 leading 
to a ratio of time scales t sys /t Tes = 16. Obviously the 
Born-Markov approximation is not valid in this parame- 
ter regime. Attenuating the coupling constant T about a 
factor of 5 suffices to reach parameters where the Born- 
Markov approximation holds. 

In the following section we apply the time- 
convolutionless projection operator technique to our 
model. This leads to a perturbative expansion for the 
master equation, which allows us to analyze the behavior 
of an atom laser in the non-Markovian parameter regime. 
As we will see in section [I D the solution of this master 



equation agrees very well with the exact solution for the 
parameters mentioned above. 



C. Time-convolutionless projection operator 
technique 

It is usually assumed that a non-Markovian dynamics 
necessarily involves a density matrix equation containing 
a time-convolution kernel. However, employing the time- 
convolutionless projection operator technique we can de- 
rive an equation of motion for the system density op- 
erator which is local in time. The following discussion 



briefly summarizes this approach. A more detailed dis- 
cussion can be found in |15|-|l7|. Especially a discussion 
of the range of validity of this approach is given in Ref. 

0- " 

The starting point is the Liouville-von Neumann equa- 
tion for the density matrix of the total system p to t 



d_ 
dt 



Ptot(t) = -ia[H int {t), p tot (t)} = aL{t)p tot {t) . (22) 



Here a denotes the strength of the coupling, and H- luij (t) 
is the interaction Hamiltonian in the interaction picture. 
If we define the projection operator V 



Vptot{t) = Tr rcs {p t ot(£)} <8> Pros = Psys(t) <S> Pr 



(23) 



it is possible to deduce an equation of motion for the 
system density operator p sys {t) 



0_ 
dt 



Vpt t(t)=K(t)V Pto t{t) 



(24) 



Here Tr res denotes the trace over the reservoir. Under 
certain conditions which are always satisfied for short 
times and in the weak coupling regime (see e.g. Ref. 17 1) 
it is possible to derive a perturbative expansion for the 
generator K(t) of the time-convolutionless master equa- 
tion in powers of the coupling strength a 

K(t)=a 2 KW(t)+a 4 K^(t)+a 6 K (6 Ht) + ... . (25) 

The superoperators K^ n \t) are given by 

t t\ tn-2 

K {n \t)= fdtlfdh... J dtn-tknit, ii...,t n _i), 


(26) 
where 



k n (t,h, . . .,i n _i) = ^(-l) 9 - 1 ^) . . . L(U) 
VL(t j )...L{t l )VL(t m )... 



(27) 



The sum in ( |27J ) is to be taken over all possible inser- 
tions of V's in between the n factors L, while keeping the 
chronological order in time between two successive inser- 
tions of V |16[ . This means t > . . . > U and tj > . . . > ti- 
The constant q is the number of inserted projection oper- 
ators P's. Assuming that all odd moments of H{ nt (t) with 
respect to p rcs vanish we have VL{t\) . . . L{t2k+i)'P = 0. 
Therefore, all terms containing odd powers of the cou- 
pling strength a disappear in equation (p5|). With the 
help of Eqs. (Eq) and (127n we obtain the superoperators 



K^ 2 \t)= / dhVL(t)L(ti)V 
Jo 



(28) 



and 



K^(t) = dt x dt 2 / dt 3 (29) 

Jo Jo Jo 

x{-pL(t)L(h)L(t2)L(t 3 )V - VL{t)L{t!)VL{t 2 )L{t 3 )V 

- VL{t)L{ti)VL{ti)L{h)V - VL(t)L(t 3 )VL(ti)L(t 2 )V} . 

We have determined also the sixth order superoperator 
K^ 6 \ Because of its length the expression for K^ 6 '(t) is 
not presented here: it is a fivefold integral of 45 terms, 
each containing 6 superoperators L(t). 

The above superoperators can be expressed in terms of 
commutators whose evaluation is simplified by invoking 
the bosonic commutation relations. For the present study 
we have determined K^ 2 \ K^\ and K^ 1 with the help 
of the computer algebra system Mathematica. This sys- 
tem enables the evaluation of the expansion coefficients 
K (n \t). 

Transforming Eqs. (3d) and (||) to the interaction pic- 
ture, the interaction Hamiltonian H m t{t) reads 

H int (t) = -ih {B(t)a){t) - B\t)a{t)} , (30) 

where 



B(t) 



dk hz(k)bke 



— iuj^t 



a(t) 



a e 



-iujot 



(31) 



The expansion parameter in the case of the atom 
laser is the coupling strength T 1 / 2 . With the as- 
sumption of an initially empty Gaussian reservoir, 
Tr res {p res i3t(i)iJ(t 1 )} = 0, we obtain a master equation 
in Lindblad form 



^Psys(t) 

1 



--S{t){<Ja,p sys (t)] 



+ 2^(*) {-^ap sys (t) - p sys (t)a t a + 2ap sys (t)a t } . (32) 

Here the Lamb shift S(t) and the decay rate ^(t) are 
given by 

S (2) (i) + S (4) (t) + S (6) (i), 



S(t) 

7 (i)= 7 (2) W+7 (4) W+7 (6) W 



(33) 
(34) 



The functions S^ n \t) and ^^(t) can be easily evaluated 
as integrals over the known functions 4>(t), ip(t) [see Eq. 
(P|)]. We find 



S<*>(t) 



ettiV(t-*i), 



dti 



dU 



and 



7 W (*) = 



(t>(t-t 2 )Tp(tl- 

•0(t-t 3 M*i- 
dti<f>(t-ti), 



dt! 



dU 



(35a) 

dt 3 [i>(t - t 2 )4>(.h - h) 

t 2 )} , (35b) 



(36a) 



dt 3 {<f>{t - t 2 )(j>{t! - t 3 ) 



+ <t>(t - t 3 )<t>(t! - t 2 ) - ij(t - t 3 )l()(t! - t 2 ) 
- iP(t - t 2 )i>(t x - t 3 )\ . (36b) 



Again we refrain from presenting the sixth order terms 
because of their length, but they are also given as simple 
integrals. 

Note that we have derived a master equation in Lind- 
blad form without making the Born-Markov approxima- 
tion. If we keep only the second orde r ter ms an d extend 
the range of integration in Eqs. ( 35a ) and ( 36a ) to infin- 
ity we get the Born-Markov master equation (|16|). Hence 
our master equation has the same simple form as equa- 
tion (Jig), but the parameter regime in which Eq. (p3) 
is valid is much larger. This will be demonstrated in 
section pD. 

From Eq. (132) one derives the equation of motion 
for the occupation number of the BEC. We are only 
interested in the normalized occupation number n(t) = 
{a' a) (t) / {a' a) (0) . One easily demonstrates that n(t) is 
given by 



n(t) = exp 



7(*i)dt] 



(37) 



For very strong coupling with T « 10 6 s~ 2 the 6th or- 
der provides reliable results only for short times. The 
rise in the atom number after the first collapse shown in 
Fig. [3| can not be reproduced by any algorithm which is 
based on a perturbative expansion because the transition 
rate diverges at this point, and the equation of motion is 
no longer analytic [p~7f. Nevertheless the short time be- 
havior is very important for the evaluation of correlation 
functions. 

Fig. Q compares the decay rates for the coupling 
strength T = 10 5 s~ 2 . As expected, the decay rate is 
valid for longer times if we include higher order terms 
in the perturbation expansion. Although the decay rate 
in 6th order does not exactly render the exact solution, 
it gives a good approximation over a wide range until 
about "fMt ~ 4. It constitutes an essential improvement 
compared to the constant Born-Markov decay rate. It is 
important to emphasize that the TCL technique is much 
simpler to perform than the solution of the exact equa- 
tion of motion. 



Hence it suffices to evaluate the integrals over the time 
dependent decay rates r )^ n '{t) in the different orders of 
the coupling strength T 1 / 2 to get the normalized occupa- 
tion number of the BEC. Provided x n '(t) is positive the 
quantity F(t) = 1 — n(t) can be interpreted as the wait- 
ing time distribution [J25 26 for a transition of a single 
atom out of the trap. 



D. Numerical results 

In our simulation we tak e th e parameters as men- 
Through variation of 



tioned at the end of section II B 



the coupling strength T we will go beyond the parame- 
ter regime in which the Born-Markov approximation is 
valid and compare the results obtained from the time- 
convolutionless projection operator technique with the 
exact solution. 

In Fig. p] the occupation number of the BEC is plotted 
for a coupling strength T = 5 x 10 4 s~ 2 . While the Born- 
Markov approximation fails to render the exact solution, 
the depicted 4th order is a very good approximation. Al- 
though it is not shown here the 6th order would be almost 
perfect. 

Fig. @ shows the results for a coupling strength T = 
10 5 s~ . Here we clearly see that for realistic param- 
eters the Born-Markov approximation fails, whereas the 
time-convolutionless projection operator technique to 6th 
order in the coupling strength reproduces the exact solu- 
tion. Our numerical results show that the Born-Markov 
approximation begins to fail at a coupling strength T w 
10 4 s -2 . This corresponds to a ratio of the time scales 
given by t sys /t YOS = 80. The TCL algorithm to 4th order 
provides reliable results until about T ss 5 x 10 4 s -2 . For 
stronger couplings to about T w 10 5 s~ 2 corresponding to 
a ratio of time scales t sys /t ICS = 8 the 6th order is needed. 



III. CONTINUOUS WAVE ATOM LASERS 

The next step after the realization of a pulsed atom 
laser was to build a continuous wave (cw) atom laser 
which was recently achieved pJq| . As in the case of the 
pulsed atom laser we model the output coupler through 
a two-photon Raman transition. Such an output coupler 
has recently been experimentally realized by Hagley et 
al. |5| . Because of the momentum transferred by the Ra- 
man transitions this is the first atom laser with a highly 
directional output. 

In addition to the pulsed atom laser discussed in the 
first part we include a pumping mechanism to describe 
a cw atom laser. In this section we extend the model 
of a pulsed atom laser to a pump mechanism already 
proposed in p-p],p7j . While in these papers the output 
coupler is treated within the Born-Markov approxima- 
tion we will use the same output coupling mechanism 
as discussed in section || to study the effects of a non- 
Markovian output coupler. The coupling to the pump 
reservoir is a fast process and may therefore be treated 
within the Born-Markov approximation. 



A. Model of a continuous wave atom laser 

We use a binary-collision atom laser scheme based on 
the near resonant dipole-dipole interaction. The atomic 
trap from which we want to study the output coupling 
consists of many atomic modes. Under certain experi- 
mental conditions it is possible to restrict ourselves to 
three modes of the trap, as discussed in pftO . Atoms 
of a thermal source - the first reservoir - are transferred 
into the pump mode of the trap. As shown in Fig. H 
these atoms enter the trap mode 1 with a rate K\N, 



They can also leave the trap from this mode with the 
rate Ki(l + N). The bosonic creation and annihilation 
operators of the first reservoir are denoted by Q\Q and 
N is the stationary occupation number of the mode 1 
if there were no transitions inside the trap. After the 
mode 1 is occupied by more than one atom, two atoms 
can interact via the near resonant dipole-dipole interac- 
tion through which one atom is transferred to the weakly 
bound mode 2 and the other to the mode of the atomic 
trap. 

The atoms in the mode 2 are out-coupled to the sec- 
ond reservoir with a rate ki . To enable evaporative cool- 
ing which is needed to achieve the necessary low tem- 
peratures for Bose-Einstein condensation the condition 
k-2 3> Ki must hold. This means that each atom which is 
transferred to the highest mode leaves the trap very fast, 
such th at th is mode is most strongly depleted. Just as in 
section [I A the atoms are out-coupled through a Raman 
transition from the ground state of the atomic trap. To 
achieve a steady occupation number of the ground mode 
the condition K\ 3> n Q must hold. This output coupling 
to the reservoir will be described taking into account 
the non-Markovian character of the dynamics, whereas 
the coupling to the other two reservoirs is treated in the 
Born-Markov approximation. This is a reasonable as- 
sumption because of the condition K2 3> ni ~3> no- 

In general we find for the binary-collision atom laser 
scheme the Hamiltonian 



H — H sys 



V + Hr + H IS + Hr 2 



(38) 



where H sys is the system Hamiltonian and V de- 
scribes the collisions between two atoms. The operators 
Hifl,Hj,i,Hi2 are responsible for the coupling to the 
three reservoirs. -ff sys and the most general form of the 
interaction Hamiltonian V can be written as 



^sys = y^?kJjOjOj , 



i=o 



V = h J^ 9j,k,l,m Ojalaia m 

j,k,l,m= 



(39) 

(40) 



In the interaction Hamiltonian V we do not consider ef- 
fects which are negligible at low densities, as the annihi- 
lation of two ground state atoms and the creation of two 
atoms in the highest modes. Such transitions are ener- 
getically unfavoured at low densities. Moreover, since 
K2 3> ki the strongly depleted mode 2 can be elimi- 
nated adiabatically ||. We introduce an effective reser- 
voir which describes the annihilation of two atoms in the 
mode 1, the creation of an atom in mode and 2, and 
the immediate loss of the atom in the mode 2. This effec- 
tive reservoir includes the terms describing collision with 
atoms in the mode 2 which were included in V . With 
these assumptions the remaining terms in V are 

V = h{ 50000 ao a o a o a o + ffim a\a\aiai 

+ ffoioi aja a{ai} . (41) 



As shown in [|9|]lC[] the first two terms in the above equa- 
tion cause a broadening of the output spectrum. In this 
paper we are only interested in the occupation number 
of the BEC, and we will not study the output spectrum. 
Because of this we can ignore these terms since V affects 
only the off-diagonal terms of the reduced density matrix 
p sys . Hence we obtain in the interaction picture 



H(t) = H It0 (t) + H Itl (t) + H oS (t), 



(42) 



where the operators describing the coupling to the three 
reservoirs are given by 



Hi ,i(t) = 
H cS (t) = 



—ih 



f 0)<4< 
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^ (QffiaU^-aiQiWe-*--' 1 



t.2 i(u -2ui)t 



, (43a) 
(43b) 



R\t)ala\e 



-a al 2 R(t)e- i( - Uo - 2un)t ) . (43c) 



The operators B,Q,R and B^,Q\R> are the annihila- 
tion and creation operators, respectively, for the reser- 
voirs shown in Fig. H. The energies of atoms in the two 
atomic trap modes and 1 are given by huo and hu>i. 
The coupling constant fi describes the strength of the 
dipole-dipole interaction. 

Summarizing, the atomic trap in which the Bose- 
Einstein condensate is built couples to three reservoirs. 
The first reservoir provides the pump mode of the trap 
with atoms from a thermal source. Through near reso- 
nant dipole-dipole interaction the atomic trap modes 
and 2 get occupied. The immediate loss from atoms in 
the second atomic trap mode to the corresponding reser- 
voir is necessary to enable evaporation cooling. Just as 
in the case of the pumped atom laser the Bose-Einstein 
condensate is coupled out from the mode through a 
two-photon Raman transition. 



B. Derivation of quantum master equations 



In this section we apply the time-convolutionless pro- 
jection operator technique to the case of the continuous 
wave atom laser. With the help of Eqs. (pq) and ( f43| ) we 
obtain to second order in the coupling strength r^ 



i:) 



- (rV) 



— Psys (t) = (£ 



+ 4 2 u } t + £ in + C coll ) Ps ys{t) . (44) 



(2) 



The superoperators describing the Lamb shift C Q and 

(2) 

the output coupling £^ u ' t to second order are defined 
through Eqs. (fig). The operators C m and C co n de- 
scribe the coupling to the thermal reservoir from which 



the atoms enter the trap mode 1, and the coupling to the 
effective reservoir which enables the evaporative cooling 
mechanism from the trap mode 2. These operators are 
treated within the Born-Markov approximation, and are 
thus independent of the expansion parameter. Hence we 
have 



4 n) a = -^ (n) (*)[«S«o.-R] 



4uU = 7 (n) W joofloj - ~aia R-^Raia 
C m R = Nk\ < a\Rai — -ai a[R — -Rai a\ > 
+ (1 + N)ki < ai Ra\ ~ -a\ a\R— -Ra\ a\ \ , 
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(45a) 
(45b) 

(45c) 

(45d) 



The func tions 7^ 2 - ) (i) and S^ 2 \t) are given by Eqs. (35a) 
and (36a). From the above master equation we obtain the 



Born-Markov approximation by replacing the time de- 
pendent functions 7 ( - 2 - > (t) and S^ 2 '(t) through the Marko- 
vian Lamb shift and the Markovian decay rate from Eq. 

©■ 
In the limit of an atom number in the Bose-Einstein 

condensate which is much greater than 1 we get the sta- 
tionary occupation number of the ground state within 
the Born-Markov approximation 



(aja ) 



Kl 



27. 



M 




(46) 



Before solving the Born-Markov master equation and the 
equation to second order perturbation theory we evalu- 
ate the 4th order term K^ from Eq. (B9T). We obtain 
after some algebra 



d_ 

dt 



— p sys (t) = ( Cq + C { J t + Cm + £coll + C oc ) Psys(t) . 

(47) 



(4) 



*(4) 



Here, the superoperators £ and £ ^t a re defined 
through Eqs. ( Eq ) and the corresp ondin g fun ction s 
7 ( 4 )(i) and S (4 >{t) are given in Eq. (§5*3) and (pb|), 



respectively. Naturally the operators £j n and £ C oii do 
not change compared to the second order expansion. In 
contrast to the second order perturbative expansion ( f44| ) 
in the coupling strength T 1 / 2 the fourth order does not 
only change the functions ^ n \t) and S^ n \i) to the ap- 
propriate order, it also adds a new superoperator C oc to 
the master equation. C oc is given by 

C oc R = r(t) ( alalRaoal — a 2 RaoaQal ) 

+r*(t) ( ala 2 Raoa\ — aoala 2 Ra\ ) 



+ - r(t) ( aoa\ a 1 RaJ ) — a aoa[ a x R ) 

+ y r *(t) ( a Rala\ 2 a 2 - Rala a\ 2 ''a\ ) . (48) 



The complex functions r(t) are defined by 



r(t) = n / dti dt 2 /(to - h) 



(49) 



The operator £ oc appears in the 4th order perturbation 

(2) 

theory due to a mixture of C ox ^ t and £ co ii- Note that 
this master equation is not in Lindblad form. Eq. ( |47| ) 
was solved by integrating the closed system of differen- 
tial equations for the diagonal elements of the reduced 
density operator. 



C. Numerical Results 

In our simulations we have chosen the parameters sim- 
ilar to the pulsed atom laser. The new parameters for the 
pumping and the coupling to the effective reservoir are 
taken to be these of Ref . |J||lO| ■ They are chosen such that 
the stationary occupation number n of the ground mode 
isn~ 100. The strength of the dipole-dipole coupling is 
il = 15 7m and the coupling constant to the pump reser- 
voir is Ki = 10 7a/. The parameter N = 20.3 describes 
the stationary occupation number of the pump mode 1 
of the atomic trap, assuming there were no transitions 
inside the trap. 

The only variable parameter in the simulation is the 
strength of the coupling T 1 ' 2 , which determines the 
strength of the coherent output coupling jm- This pa- 
rameter is varied in our simulations from a value where 
we expect the Born-Markov approximation to hold, to 
parameters where the time-convolutionless projection op- 
erator technique for the pulsed atom laser was valid in 
the first part of this paper. 

Fig. g shows the occupation number of the Bose- 
Einstein condensate for a value T = 5 x 10 4 s~ 2 . For this 
coupling strength the TCL algorithm agreed very well 
with the exact solution in the case of the pulsed atom 
laser. In the case of the cw atom laser we see strong 
oscillations in the occupation number of the BEC. Note 
also, that in 4th order the mean occupation number in- 
creased compared to the stationary occupation number 
resulting from the Born-Markov approximation. 



If we reduce the coupling strength T to T = 1 x 10 4 s -2 , 
then there are only very weak oscillations around the sta- 
tionary value obtained from the Born-Markov approxi- 
mation and the fourth order perturbation theory agrees 
well with the second order expansion. In the case of the 
pulsed atom laser the Born-Markov approximation also 
failed for coupling strengths T > 1 x 10 4 s -2 . 

As expected this suggests that for the cw atom laser the 
parameter regime in which the TCL-algorithm is valid 



coincides with that of the pulsed atom laser. Thus the 
parameters in Fig. y are chosen such that the TCL algo- 
rithm to 4th order should provide reliable results. 

The depicted oscillations in the atom number can be 
interpreted as an interference effect between two favored 
transition modes. The physical origin of this interference 
phenomena can be seen already from the second order 
term 7^ (t) which takes the form 



7 (2) (i)= / ds(f>(s) 
Jo 



= 2 



dm J(uj) 



sin[(wo — w)i] 

LUq — UJ 



(50) 



Here the first factor of the integrand, namely the spec- 
tral density J(lu), has a uj^ 1 ^ 2 singularity at w = 
[see Eq. (0)], whereas the second factor is concentrated 
around the ground state frequency lo = loq. Thus, the 
main contribution to this integral stems from transitions 
with frequencies near uj = and lo — loq. The corre- 
sponding transition amplitudes interfere and lead to the 
observed oscillations in the decay rate and therefore to 
the oscillations in the occupation number. 

In a Markovian system, where the spectral density of 
the coupling strength J(ui) is bounded, the second factor 
of the integral (pQ) quickly approaches a S- function in the 
limit t — > 00, and the oscillations decay on a time scale 
i res . This is known as Fermis golden rule. 

In contrast, in the case of an atom laser, the contribu- 
tion near u> = is also important for long times. This 
can be seen by the asymptotic behavior of 7 < - 2 - ) (t) which 
can be evaluated explicitly from Eq. ( |50| ) and is given 
through 



7 (2) W=7M 



2F 



A/awp 



cos(ujQt + 7r/4) 



(51) 



According to this relation the decay rate 7^ 2 - ) (£) ap- 
proaches the Markovian rate 7m very slowly as t~ x / 2 . 
It is important to note that this behavior is due to the 
uj^ 1 ' 2 singularity of J(lo) at uo = which, in turn, is a di- 
rect consequence of the dispersion relation (0) of massive 
particles. Thus we observe that non-Markovian effects 
decay only slowly in time with an algebraic behavior. 
This clearly demonstrates that the golden rule limit 7j\/ 
is relevant only for extremely long times and shows the 
importance of a non-Markovian treatment of the atom 
laser. 

From Eq. ( J50| ) one clearly sees that the oscilla- 
tions will also appear if we abandon the assumption 
fco ~ 0. The function J(u>) will then be peaked around 
ljl = hk 2 /(2M) but still has a u -1 / 2 singularity at 
U) = 0. Hence a significant contribution to the inte- 
gral in Eq. ( p0| ) comes from u> — and, therefore, the 
time dependence of 7^ 2 '(i) again shows an oscillatory be- 
havior. These oscillations are also present in the case of 
the damped Jaynes-Cummings model with detuning (see, 
e.g., Ref. |17]]), but there they are damped exponentially. 



IV. SUMMARY 



We have demonstrated that in the case of a pulsed 
atom laser with a Raman output coupler in a realistic 
parameter regime in which the Born-Markov approxima- 
tion fails, the time-convolutionless projection operator 
technique is able to reproduce the exact solution. 

If we leave the intermediate coupling regime and take 
the 6th order terms in the output coupling strength we 
are able to render the exact solution for short times to 
about t = 4/7m in the strong coupling regime. This 
short time behavior is important for the evaluation of 
correlation functions. The long time behavior after the 
first collapse of the atom number can not be reproduced 
with the TCL algorithm, because the decay rate diverges 
at such a point. 

In the case of a cw atom laser we have chosen the cou- 
pling strength such that the perturbative expansion for 
the pulsed atom laser agreed with the exact solution for 
the considered times. We found strong oscillations of the 
occupation number of the BEC and the mean atom num- 
ber increased. Our results clearly suggest that these os- 
cillations survive for very long times because of the slow 
algebraic decay of the relaxation rate. 
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FIG. 1. Normalized occupa- 

tion number n(t) = (a t a)(t)/{a t a)(0) of the BEC in the 
atomic trap. Plotted are Born-Markov approximation, ex- 
act solution and perturbation expansion to 4th order in the 
coupling strength F = 5 x 10 4 s~ 1 . This leads to a ratio of 
time scales t EyB /t TCE — 16. 



FIG. 2. The normalized occupa- 

tion number n(t) — (a t a)(t)/(a t a)(0) of the BEC: exact 
results and the results from the Born-Markov approxima- 
tion, resp., from 4th and 6th order perturbation theory. 
The parameter F = 10 5 s -1 is chosen such that the ratio of 
time scales becomes t ras /t aya — 8. 



exact solution 

TCL6 

Born-Markov approximation 





FIG. 3. Normalized occupation num- 

ber n(t) = (a t a)(i)/(a t a)(0) of the BEC: Born-Markov 
approximation, exact solution and perturbation expansion 
to 6th order in the coupling strength F = 10 6 s -1 . This 
leads to a ratio of time scales t las /t sys = 0.8. 



FIG. 4. Comparison of decay rates. Plotted are 
Born-Markov, exact decay rates and the decay rates in 
4th and 6th order perturbation theory. The value of the 
coupling strength is T — 5 x 10 4 s _1 . 
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FIG. 5. Our simplified model of a cw atom laser consists 
of a pump mode 1, the ground state mode which is oc- 
cupied from the Bose-Einstein condensate and the weakly 
bound mode 2 from which the atoms are strongly coupled 
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FIG. 6. Occupation number n(t) — (a/a)(t) of the BEC 
in an atomic trap in the case of a continuous wave atom 
laser. Plotted are the Born-Markov approximation and the 
results of a non-Markovian perturbative expansion of the 
output coupling to second and fourth order in the coupling 



